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Abstract 

In the presented paper we raised the question whether initial cosmological singularity can be 
proved by cosmological tests. The classical general relativity theory predicts the existence of 
singularity in the past if only some energy conditions are satisfied. On the other hand the latest 
quantum gravity applications to cosmology suggest the possibility of avoiding the singularity and 
replacing it with a bounce. Bounce is the moment in the evolution of the Universe when the 
Universe's size has minimum. Therefore the existence of observationally detected bounce in past 
of Universe could indicate the validity of the loop quantum gravity hypothesis and nonexistence 
of initial singularity which is present in the classical ACDM. We investigated the bouncing model 
described by the generalized Friedmann-Robertson- Walker (FRW) equation in the context of the 
observations of the currently accelerating universe. The distant type la supernovae data are used 
to constraint on bouncing evolutional scenario where square of the Hubble function H 2 is given 
by formulae H 2 = Hj 2 [Q m fi(l + z) m — fi n ,o(l + z) n ], where fi m ,o, Q n ,o > are density parameters 
and n > m > 0. In this paper are showed that the on the base of the SNIa data standard bouncing 
models can be ruled out on the 4cr confidence level. After adding the cosmological constant to 
the standard bouncing model (the extended bouncing model) we obtained as the best-fit that 
the parameter fi raj o is equal zero which means that the SNIa data do not support the bouncing 
term in the model. The bounce term is statistically insignificant on the present epoch. We also 
demonstrated that BBN offers the possibility of obtaining stringent constraints of the extra term 
f2 n ,o- The other observational test methods like CMB and the age of oldest objects in the Universe 
are also used. We use as well the Akaike informative criterion to select a model which fits data the 
best and we concluded that bouncing term should be ruled out by Occam's razor, which makes 
the big bang scenario more favorable then the bouncing scenario. 

PACS numbers: 98.80.Bp, 98.80.Cq, ll.25.-w 



2 



I. INTRODUCTION 



We are living in an age of high precision cosmology which offers a possibility of testing 
exotic physics, which is obvious for the early Universe In this context the most im- 
portant are BBN constraints because the present Universe opens only small windows on 
the exotic physics. The main aim of this paper is to discuss whether the initial singularity 
can be checked against the astronomical observations. The question of singularity cannot 
be answered directly, therefore we use two prototype models based on the classical and 
quantum gravity theory. The first is the ACDM which is a concordance model describing 
the evolution of the Universe from the initial singularity (the big bang) driven by the cold 
dark matter and the cosmological constant (dark energy). The second is a bouncing model 
which appears in the context of quantum cosmology and characterized by the lack of initial 
singularity. During its evolution, the expansion phase is proceeded by the contraction phase 
at the bounce where the scale factor assumes the minimum nonzero value. 

We use some tests to discriminate between these two alternative models. One of the most 
important tests applies the SNIa data to fit the cosmological models. Recent measurements 
of type la supernovae observations suggest that the universe is presently accelerating 
A dark energy component has usually been proposed as a source of acceleration mechanism 
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Many theoretical propositions have been suggested about these components. However, 
the different effects arising from quantum fluctuation, spinning fluid, etc. can also mimic 
dynamically the role the dark energv which drives acceleration through an additional term 



in the Friedmann equation 
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] . Some of them give rise to the 



bounce. In many cases they prevailed in the very early epoch but are very small in the 
present epoch. Therefore it is very difficult to detect the existence of this component in the 
present and those of the relatively close past (after CMB) observations of SNIa. 

In the present work we investigate observational constraints on the evolutionary scenario 
of the standard bouncing cosmological models defined as a class of models for which the 
Hubble function H and the scale factor a are related by the formula 

h 2 = H 2 (n mfiX - m - n nfiX ~ n ) (i) 

where n > m > and x = ^ where the index zero denotes the quantities evaluated in the 
present epoch, the parameter f2 n is called the bouncing term; and the density parameters 
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satisfy the constraint relation 

^m,0 — ^n,0 = 1- 

While focus mainly on the constraints coming from SN la data and WMAP observations, 
the complementary constraints coming from BBN and the age |l5( of the oldest high-redshift 
objects are also considered. We use the maximum likelihood method to estimate the model 
parameters m, n and Q n ,o- Similarly we analyze the models with the additional parameter — 
the cosmo logical constant. It is called the generalized bouncing model. 

The proposition of the bounce type evolution of the early universe seems to be very 
attractive not only from the point of view of the quantum description of the early Universe 
because the expansion of the universe is accelerated automatically due to the presence of 
the bouncing term. 

The standard bouncing scenario predicts the acceleration around the bounce with a tran- 
sition to the deceleration epoch. The cosmological constant brings this deceleration epoch 
to the end and a new acceleration epoch begins. 

Therefore, these models can be proposed as the models of our Universe, because they 
include the epoch of acceleration. However we show that the influence of the bouncing term 
is insignificant in the present epoch. Therefore, the data from the present epoch, such as the 
SNIa data, have not power to rsider the model with bouncing term statistically significant. 
So the ACDM model with big-bang scenario is strongly favored by data over the model with 
the bounce. 

By the application of standard Akaike criterion of the model selection we can choose 
the ACDM model over the generalized bouncing model. We conclude that the data fail to 
support the existence of the bouncing term. The bouncing term in the present epoch is 
insignificant and it is not possible to detect its influence by the use the latest SNIa data. 

This fact justifies certain scepticism about the existence of the SNIa window on exotic 
physics in the current epoch. However we cannot rule out other models by testing them 
against the current data. It is also possible to investigate the differences in the predictions 
of these models for some earlier epoch. 

For example the BBN epoch is a well tested area of cosmology. From this analysis 
we gather that the extra term Q n ^x~ n causing the bounce should be constrained to be 
sufficiently small during nucleosynthesis. 

The organization of the text is the following. In section 2 the evolutionary scenario of 
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bounce FRW cosmologies is investigated by the use of dynamic system methods. We show 
that they are structurally unstable due to the presence of centers in the phase portraits. In 
section 3 we discuss the constraints from SNIa data on the standard bouncing models. In 
section 4 we extend the bouncing models by introducing the cosmological constant and then 
we study how these models fit the current supernovae and WMAP data. In section 5 we 
formulate conclusions. 



II. THE BOUNCING MODELS: BASIC EQUATIONS 

The idea of bounce in FRW cosmologies appeared in Tolman's monograph devoted to 
cosmology Q]. This idea was strictly connected with oscillating modei 
present oscillating models play an important role in the brane cosmology 
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universe undergoing a bounce instead of the big-bang is also an appealing idea in the context 



of quantum cosmology 



22| . The attractiveness of bouncing models comes from the fact that 



;hey have no horizon problem and they explain quantum ori gin of structures in the Universe 
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. Molina-Paris and Visser and later Tippett j2y, |2jj characterized the bouncing 
models by the minimal condition under which the present universe arises from a bounce from 
the previous collapse phase (the Tolman wormhole is different name for denoting such a type 
of evolution). The violation of strong energy condition (SEC) is in general a necessary (but 
not sufficient) condition for bounce to appear. For closed models it is sufficient condition 
and none of other energy condition need to be violated (like null energy condition (NEC): 
p + p > 0, week energy condition (WEC): p > and p + p > 0, dominant energy condition 
(DEC): p > and p ±p > energy conditions can be satisfied). 

We can find necessary and sufficient conditions for an evolutional path with a bounce by 
analyzing dynamics on the phase plane (a, d) , where a is the scale factor and dot denotes 
differentiation with respect to cosmological time. We understand the bounce as in 

BEL 

namely there must be some moment say t = tbounce in evolution of the universe at which 
the size of the universe has a minimum, db OU nce — and a > 0. This weak inequality a > 
is enough for giving domains in the phase space occupied by trajectories with the bounce. 
Let us consider the dynamics of the FRW cosmological models filled by perfect fluid with 
energy density p and pressure p parameterized by the equation of state in the general form 



p = w(a)p. 



(2) 



The basic dynamical system constitutes two equations 

! = -£(, + %,) (3) 
P=-3H(p + p) (4) 

Equation (jSJ) is the Rauchadhuri equation while equation (jlj) is the conservation condition. 
If the equation of state is postulated in the form (J2J then from (j3J) we obtain 



a w(a) 



p = p(a) = p a exp ^—3 J da J (5) 

It is interesting that dynamics of the model under consideration can be represented in the 
analogous form to the Newtonian equation of motion 

a = -da- (6) 

2 

where V = plays the role of potential function for the FRW system. Therefore different 
cosmo logical models are in a unique way characterized by the potential function V = V(a) 
and we can write down the Hamiltonian for the fictitious particle-universe moving in the 
one-dimensional potential as 

2 

H = ^ + V(a), p a = a (7) 
It is useful to represent eq. (JHJ) in the form of dynamical system 

x = y, y = ~^, (8) 
where we denote x = a, y = a and system (jHJ) has the first integral in the form 

y + V(x) = -~ (9) 



where k is the curvature index. 

The critical points of the system (jHJ) if exist are: yo = 0, (f^)x = 0, i.e. they are always 
static critical points located on x-axis. The form of first integral Q defines the algebraic 
curves in the phase plane (a, a) on which lies solutions of the system. This solutions are 
in two types: regular is represented by trajectories or singular is represented by singular 
solutions for which the right-hand side of (jHJ) are null (or V{xq) = — | for nonflat models). 
Note that the bouncing points are intersections points of trajectories situated in the region 
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of the configuration space in which < 0, i.e. V(a) is a decreasing function of a or has 
extrema. It is well known that the systems in the form (JHJ) have only critical points of two 
admissible types: centres if (^pr)x > or saddles in opposite case if (^■) XQ < 0. Therefore 
all trajectories with bounce intersect an x-axis and then they are situated on the right side 
from the critical point at which a = and a > 0. The critical points are represented by 
points as well as by separatrices of the saddle point. In others words bouncing trajectories 
are represented by such trajectories in the phase plane which are passing through the x-axis 
in such a direction that they always belong to the accelerating region (in the neighborhood 
of bounce). Of course it is only possible if the SEC is violated. 

Let us consider some prototype of bouncing models given by the Friedmann first integral 
in the form 

h 2 = 4 - 4> (io) 

a m a n 

where A, B are positive constants and n > m, H = (In a)' is the Hubble function and a dot 
denotes differentiation with respect to cosmological time t. 
It is convenient to rewrite (|T0|) to the new form 

H 2 = H 2 (Q mfi x- m - n n , x- n ), (11) 

where fi m0 , Q n ,o are density parameters for noninteracting fluids which give some con- 
tributions to right-hand sides of eq. (fTUj) . We define density parameters f2 mj0 = ™ ' 



'0 

^no — ^%tt~i where an index "0" means that corresponding quantities are evaluated at the 
present epoch, x = is the scale factor expressed in the units of its present value ao- 

After differentiation the both sides of (jll)) with respect to the reparameterized time 
variable r: t — > r, \H Q \dt = dr we obtain 

- = \ (a m , (2 - m) x~ m + n nfi (n - 2) x~ n ) (12) 

If we consider the generalization of the bouncing models with the cosmological constant 
then in both equations (fTTj) and (fT2j) the parameter fi^o should be added to their right-hand 
sides. 

Note that the bouncing models can be treated as the standard FRW models with two 
noninteracting fluids with energy density and pressure in the form 

P = p m + Pn = 3H$tt mfi x~ m - 3H$tt nfi x~ n 
( m\ ( n\ 

P = y-l + -j) Pm + [-I + g J Pn Pm > 0, p n < 0. 
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The curvature term as well as cosmological constant term can be obtained in an analogous 
way by putting m = 2 or m = 0, respectively. 

If we postulate that the present universe is accelerating, i.e. 'x > at x = 1 then in the 
general case with the cosmological constant we obtain the following condition 

fi m , (2 - m) + n nfi (n - 2) + 2ft A , > (13) 

Because relation (fTHJl is validate any time, the substitution H = H Q and x = 1 to (fT3j) gives 
constraint 

^m,0 — ^n,0 + ^A,0 = 1 (14) 

Let us now concentrate on the standard bouncing models (SB) without the cosmological 
term. Then from (J13j) including constraints (|14|) we obtain the sufficient condition for ac- 
celeration at present 

Qm,o( n — to) > 2 — n, (15) 

where for the case of n = 2, k = 1 fi m .o > is only required for present acceleration. 

If only f2„ o is larger then Q m ,o the bouncing universe is presently accelerating for any m, 
n parameters. It is worthy to mention that condition (|T3jl is minimal qualitative information 
about acceleration and the rate of this acceleration is required for explanation SNIa data. 

From the definition !) 11)1 one can obtain the domain admissible for motion of the bouncing 
models 

D = {x: x > Xb where x^ = I w '° ) }. (16) 



From 1)11)1 the potential function V(x) in the particle-like description can be determined 

Vix) = = ~\n eff {x)x* (17) 

where effective density parameter Q e ff — ^m.o x ~ m ~~ ^n,o x ~ n - The acceleration region in 
the phase plane can be determined in term of potential function, namely if 

then universe is accelerating. 

From !)17p we obtain result that at the bounce moment 



fdV\ 1 



\dx 2" hU \n 

\ / x b \ 



* = -(^r) =^o(S^) (n-m) (19) 



'n,0 



which indicate that bouncing models defined by equation (JTTJl at the bounce are in accel- 
erating phase for any ranges of model parameter m, n, Q m ,o, ^n,o- Because a& is larger then 



a : (^)ao = the bouncing universe stay still in the accelerating region. 

From eq. we obtain that the universe start to accelerate at the point x = xo such 
that 



i 

m — n 



^ = n (m L (2Q) 



Q nfi (n - 2 

where a positive value of (to — 2)(n — 2) is required. Note that in any case 

x < x fe (21) 

i.e., the start of acceleration proceeds the bounce. The value of xq determine the location of 
the critical point of the dynamical system x = y, y = — |p on x-axis. The sign of the second 
derivative of potential function determines the type of critical points (centre or saddle) 
because the eigenvalues of the linearization matrix of the system satisfy the characteristic 
equation 



From (JT7j) we obtain 

'd 2 V\ , , I 



A2 + I % ' {Xo) = °" (22) 



2 v x ) = --fin,oXo" [(2 - m){\ - to) - (2 - n)(l - n)] . (23) 

Therefore if (to, n) belong to the interval (3/2, oo) then we have centres, while if (to, n) 
belong to interval (— oo,3/2) we obtain saddles. 

Let us concentrate, for example, let us concentrate on the case of m = 3 (dust matter). 
Then 

(x ) = ±n n , x^n(n - 3). (24) 

Hence, if n > 3 we obtain (^^j (%o) positive which corresponds to the centres in the phase 
plane. The presence of centres on the phase portraits means that all bouncing models are 
oscillating. Let us note, however that the corresponding systems are structurally unstable 
because of the presence of nonhyperbolic critical points on the phase portraits. Physically 
it means that the small perturbation of right-hand sides of the system under consideration 
disturbs a qualitative structure of the orbits (i.e. a phase portrait). On Fig. [T] and 121 
the phase portraits and diagrams of the potential function are presented for the standard 
and generalized bouncing model. Fig. [T] describes all special cases listed in Table HI The 



TABLE I: Some special cases of bouncing models 



model 


dynamical equations (first integral) 


FRW model dust filled 
universe with global rotation [14] 
or brane models with dark radiation [40] 


x = yx 3 

V = h ( — ft™ nX~ 2 + 2lf2, , n\x~ S ) X 3 

i/ 2 V III , U 1 "J u "(jj jU | **-* / 

= \ (-Qrn fi X~ l - \Quj,o\ X ~ 2 + fifc.o) X 6 


FRW dust filled 
universe with spinning fluid [52] 
or a class of MAG models [301 


x = yx 5 

ii = 1 f — Q n x -2 + 40, nlaf 5 ') x 5 
^ = i (-il mi0 ^ _1 - |O Si0 |x" 4 + n fc)0 ) x 10 


Stephani models 
filled by perfect fluid 
P = IP l"13] 


x = y 

y = i (-O 7)0 (l + 3 7 )x- 2 - 3 ^ + Sn^x 5 - 1 ) 

£ = \ (^ox- 1 - 3 ^ + n kfi x s ) 



classically forbidden region for a < ao is shaded. The evolution of the model is represented 
in the configuration space by a Hamiltonian level 

H = E = -Qk,o- 

The trajectory of the flat model separates the regions occupied by both closed and open 
models. The decreasing of the potential function with respect to the scale factor deter- 
mines the domain of phase space occupied by accelerating trajectories. The bounce is the 
intersection point of trajectory with the axis a. Note that around the bounce we have ac- 
celeration. On the phase plane of Hamiltonian dynamical systems only centres and saddle 
points are admissible. The centres are structurally unstable while saddles are structurally 
stable. Because the centre appears in the phase portraits of standard and generalized bounc- 
ing models, both models are structurally unstable. The critical points represent the static 
universes. Generalized bouncing model has two disjoint acceleration regions. The first is 
due to bouncing term while the second is forced by the cosmological constant term. 

The full knowledge of the dynamics required its analysis at infinity, i.e. at the circle at 
infinity x 2 + y 2 = oo. The standard procedure is to use projective maps on the plane and 
then analyze the system in a standard way. One can find the critical points at infinity as 
an a intersections of trajectory of a flat model with circle at infinity, i.e. {{x,y) : Qk,o — 0} 
and {z = l/x,u = y/x, and x = oo,v = 1/y, w = x/y, y = oo} - the trajectory of the 

2 

flat model \ = —V(x) with a circle at infinity. 
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bouncing point 



open 
models 




FIG. 1: The phase portrait and the diagram of the potential function for BM model (all case from 
Table [J). The minimum of the potential function corresponds to a centre on the phase plane. The 
acceleration region is located on the right from the ftmin 
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bouncing point 




FIG. 2: The phase portrait and the diagram of the potential function for ABCDM model. The 
minimum (maximum) of the potential function corresponds to a centre (saddle) on the phase 
plane. The system is structurally unstable because of the presence of nonhyperbolic critical point 
(a centre). 
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Some special cases of the bouncing systems contain Table |l] Because we are dealing with 
autonomous dynamical systems their phase portrait is always given modulo diffeomorphism 
or equivalent modulo any time reparameterization following the rule: r — > rj: dr = f(x)drj, 
where f(x) is diffeomorphism, x is point which belong to phase plane. 

For analysis of bouncing models in term of dynamical systems, it is useful to reparam- 
eterized original time variable in order to obtain nondegenerate critical points at infinity. 
Then we obtain -jfe and \ = |(f2 mj0 a; 2 ~ m+/3 — Q n fiX 2 ~ n+f3 ) is now representing the trajec- 
tory of the flat model, f3 should be chosen in the suitable way to regularize critical points. 
Let (3 = n — 2 then as x goes to infinity then y also goes to infinity. Only the sign of the 
parameter m (if m < 0) decides whether the future of the system is the type of a big rip 
singularity. If m = then the case of cosmological constant can be recovered. 

It is convenient to regularize the system by multiplication both sides of the system x 3 in 
the first case and x 4 in the second one respectively. It is equivalent to reparameterized time 
variable following rule r — > rj : % — dr], where f3 — 3 and (5) for the system from Table|l] For 
both systems from the table we have an additional term in the generalized FRW equation. 
In the first case the effects of global rotation produce contribution corresponding to the 
negative energy scaling like radiation. The same contribution appears in the brane models 
on the charged brane. It is known as the dark radiation [10|. Please note, that analogous 
term appeared if we include the Casimir effect coming for example from quantum effects of 
massless scalar fields 22, 28, 291. In the second case (Table HI) it is the model with spinning 

nn 

dust fluid. It can be also recovered as a class of MAG models j.7, 30]. 

In both cases we can find centre at finite domain and periodic orbits. At infinity we 
have unstable and stable nodes at x — +oc,y = =)=oo. The trajectory of the flat model 
separates the regions occupied by closed and open models. All models have bounce but 
some from them are oscillating models without the initial and final singularities. For our 
future investigations of observational constraints on bouncing models it is convenient to 
derive crucial formulae for H(z) where z is redshift z: 1 + z — x~ l . We obtain from (fTUj) 
that H 2 = Hq [Q TO) q(1 + z) m — Q nt o(l + z) n }. It is useful to represent it in the corresponding 
bouncing parameters. 

For this aim we find x\, corresponding to the bounce and value of redshift which identify 

i i 

this moment during the evolution: Xf, = J " ™, z& = — 1 + (^jpf 1 " ™ Finally we 

obtain independent model parameters characterizing its role in evolution (modulo present 
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value of Ho), namely 

H = H 




If Q m ,o is fixed, for example from independent galaxy observations then the evolutionary 
scenario is parameterized by single n parameter 




H = H ^ (1 + zf 2 J 1 + 1 - (1 + z)\ (26) 



^3,0 

where we put Q mt o = f2 3i0 , i.e. dust filled universe. 

In the case of generalized bouncing models the potential function takes the following form 

v(x) = hi nfl x 2 - n - Ul mfl x 2 ~ m - hi kfiX 2 . 

It means that if only n, m > then we obtain the de Sitter solution as a global attractor 
in the future. In the opposite case if m > the big-rip singularities are generic future of 
the model. Note that in the class of generalized bouncing models only trajectories around 
point (xq, 0) represent oscillating models without a singularity and there is admissible large 
class of closed, open and flat models which evolve to infinity. 

It is interesting that the characteristic bounce can be defined in terms of geometry of 
potential function only. By bouncing cosmology we can understand all cosmological models 
for which the potential function has at some point a minimum. 

III. BOUNCING MODEL AND DISTANT SUPERNOVAE OBSERVATIONS. 

In this section we confront the bouncing cosmological models with observations of distant 
SNIa. These observations in the framework of the FRW model indicate that present accel- 
eration of our Universe is due to an unknown form of matter with negative pressure called 



dark energy 3]. Apart from the cosmological constant t 
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rere are also other candidates for 
32, 3^. We use the SNIa data to 



dark energy which were tested from SNIa observations 
test the acceleration in the bouncing models. Moreover these models are attractive because 
they have no horizon and initial singularity, and they yield an explanation of structures 
which originated in the quantum epoch j^]. 

We consider the flat FRW model since there is a very strong evidence that the Universe 
is flat in the light of recent WMAP data js^J. We confront the two "bouncing" models (with 
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and without extra A fluid) with SNIa data. For this purpose we calculate the luminosity 
distance in a standard way 

d 'M = i l+ *£m (27) 

To proceed with fitting models to SNIa data we need the magnitude-redshift relation 

m(z,M, Q m ,o, ^a,o, n, m) - M = M + 5 log 10 D L (z, Q mfl , ^a,o, m, n) (28) 
where M being the absolute magnitude of SNIa and 

D L (z, n mi0 , Q A)0 , m, n) = H d L (z, H , tt mfi , Q A>0 , m, n) (29) 

is the luminosity distance with H factored out, so that marginalization over the parameter 
M 

M = -5 log 10 H + 25 (30) 

reads actually marginalization over H . 

The parameter A4 is actually determined from the low-redshift part of the Hubble dia- 
gram which should be linear in all realistic cosmologies. It lead to value of Hq ~ 65 km/s 
Mpc QEIQ, i.e., A4 ~ 15.955. In further analysis we estimate the models with this value 
of M. and without any prior assumption on H Q . 

Then we can obtain the best fit model minimizing the function x 2 

f,,theor ,,obs\2 

x , = j2 (jh -ft ) (31) 

i 1 

where the sum is over the SNIa sample and crj denote the (full) statistical error of magnitude 
determination and /z; = mi — Mj. 

Because the best-fit values alone are not sufficient, the statistical analysis is supplemented 
with the confidence levels for the parameters. We performed the estimation of model pa- 
rameters using the minimization procedure, based on the maximum likelihood method. We 
assume that supernovae measurements came with uncorrelated Gaussian errors and the 
likelihood function C could be determined from the chi-square statistic £ cx expf— y 2 /2) |2j. 

The first published large samples of SNIa appeared at the end of the 90s (2, |3j. Later 



other data sets have been made either by correcting errors or 
The latest compilation of SNIa was prepared by Riess et al. 



jy adding new supernovae. 
3511 and became de facto a 



standard data set. It should be noted that this compilation encloses the largest number of 
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TABLE II: Results of the statistical analysis of the bouncing model without dust (BM) and bounc- 
ing cold dark matter model (BCDM) obtained for SNIa data from the best fit with minimum x 2 
(denoted as BF) and from the likelihood method (denoted as L). The case of a fixed value of the 
parameter A4 is denoted as F. If in BF method we obtain f2 nj o = than n could be taken arbitrary 
(marked as ^4). 



model 


^m,0 


m Q n fi n M. x 2 


method 


BM model 




1.00 1.4 0.00 A 15.975 181.6 


BF 






1.00 1.5 0.00 1.7 15.975 


L 






1.54 1.4 0.54 1.5 F15.955 182.3 


BF 






1.00 1.4 0.00 1.6 F15.955 


L 


BCDM model 




1.86 


- 0.86 3.7 16.085 217.4 


BF 


(dust matter m = 


3) 


1.86 


- 0.86 3.7 16.095 


L 






1.86 


- 0.86 3.7 F15.955 273.7 


BF 






1.86 


- 0.86 3.7 F15.955 


L 



high-redshift z > 1.25 objects in compare to older compilations. From this compilation we 
take the "Silver" sample which contains all 186 SNIa, and the restricted "Gold" sample of 
157 SNIa (with higher quality of the spectroscopic and photometric records). 

In order to test a cosmological model we calculate the best fit with minimum \ 2 as weu 
as estimate the model parameters using the maximum likelihood method 2\. For both 
statistical methods we take the parameters m and n in the interval [0, 10], n > m. We test 
separately the models with and without the cosmological constant term. We also assume 
priors about fl™n and we estimate it or take f2 mi o = 0.3 (baryonic plus dark matter in 
galactic halos) [36]. 

The results of two fitting procedures performed on the "Gold" sample for the cosmological 
bouncing models with different prior assumptions are presented in (Table ITT1 and IHlj) . These 
tables refer both to the \ 2 (best fit) and results from marginalized probability of density 
functions. 

At first we analyzed bouncing model without any priors for m parameter (BM). We obtain 
the value x 2 = 181.6 what means that this model is acceptable on the 2a level with degree 
of freedom df = 153. However, the estimated value of m = 1.4 in the model is unrealistic 
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TABLE III: The results of statistical analysis of BCDM models (m = 3) obtained for SNIa data 
from the best fit with minimum x 2 (denoted as BF) and from the likelihood method (denoted as 
L). The case of a fixed value of A4 is denoted as F. 



model 


n fin fl A4 V 2 


method 


BCDM model 


1.50 0.50 16.105 226.6 


BF 


with. 7i = 4 


1 50 50 16 095 — 


L 




1.50 0.50 F15.955 296.4 


BF 




1.50 0.50 F15.955 


L 


BCDM model 


1.03 0.03 16.175 291.2 


BF 


with n = 6 


1.03 0.03 16.175 


L 




1.03 0.03 F15.955 443.4 


BF 




1.03 0.03 F15.955 


L 



because the dust matter is present in the universe (m = 3). With the prior m = 3 we obtain 
X 2 = 217.4 with the value of the parameter n = 3.7. For the more realistic model with 
m = 3 and n = 4 (because of the presence of radiation matter in the Universe) (Table IIII|) 
we obtain \ 2 = 226.6. While the bouncing model with dust (BCDM) is better fitted than 
the Einstein-De Sitter model it is rejected at least on the 4<r level. With priors M. ~ 15.955 
the model is rejected on the 8<r level. 

In Fig. El we present of residuals plots of the m-z relation for considered models with 
respect to the Einstein-de Sitter (CDM) model. Apart the CDM model (the zero line) 
the three models ACDM, BM and BCDM are shown. The diagrams for bouncing models 
intersect the ACDM diagram in such a way that the supernovae on intermediate distances 
are brighter then expected in the ACDM model, while very high redshift supernovae should 
be fainter then they are expected in the ACDM model. Note that this effects are more 
stronger for the BCDM model than for the BM model. 

Similarly we analyze the bouncing models with the additional parameter — the cosmolog- 
ical constant. We fixed value of m = 3 (dust matter) and it is called the extended bouncing 
model (ABCDM). This model with df = 153 is statistically admissible on the 2a level (Ta- 
ble E}, but we obtain O ni o = (no bounce term) as a most probable value. This result is 
similar also for models with fixed values n = 4 and n = 6 (Table IVj). as well as it is inde- 
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FIG. 3: Residuals (in mag) between the Einstein-de Sitter model and the Einstein-de Sitter itself 
(zero line), the ACDM flat model (upper curve) the best-fitted BM model, (upper-middle curve) 
and the best-fitted BCDM model with m = 3 (lower-middle curve) (with assumed M=15.955). 



pendent from the assumption on fi mi o- In this way ABCDM reduces to "classical" ACDM 
model. 

The confidence levels in the (Q n fi,n) plane are presented in Fig. 0] In order to complete 
the picture we have also derived one- dimensional probability distribution functions (PDF) 
for f2 n (Fig. |3J) and n (Fig. lUJ) obtained from the joint marginalization over remaining 
model parameters. The maximum value of such a PDF informs us about the most probable 
value of fi n ,o> supported by supernovae data within the extended bouncing dust model. 

From the PDFs the most probable value of Q Uj o is also equal 0, however non-zero value 
of Q n fi cannot be excluded. In this way, it is crucial to determine which combination of 
parameters give the preferred fit to data. This is the statistical problem of model selection 
371 ] . The problem is usually the elimination of parameters which play insufficient role in 
improving the fit data available. Important role in this area plays especially the Akaike 
information criterion (AIC) [2^. This criterion is defined as 



AIC = -21n£ + 2fc 



(32) 
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TABLE IV: Results of the statistical analysis of the extended bouncing models (m = 3), obtained 
for SNIa data from the best fit with minimum x 2 (denoted as BF) and from the likelihood method 
(denoted as L). The case of a fixed value of parameter Q m o is denoted as F. If in BF method we 
obtain f2 n) o = than n could be taken arbitrary (marked as A). 



ttioHpI 




^A,0 


M x 2 


mptTioH 


ABCDM model 


SI D DO A 

U.(J1 \J • \J\J -ii. . 


0.69 15.955 175.9 


BF 




0.34 0.00 3.0 


0.67 15.965 


L 




F0.30 0.00 A 


0.70 15.955 175.9 


BF 




F0.30 0.00 3.0 


0.70 15.945 


L 


ABCDM model 


0.31 0.00 A 


0.69 


- 175.9 


BF 


with M = 15.955 


0.31 0.00 3.0 


0.68 




L 




F0.30 0.00 A 


0.70 


- 175.9 


BF 




F0.30 0.00 3.0 


0.70 




L 



where £ is the maximum likelihood and k is the number of the parameter of the model. 
The best model is the model which minimizes the AIC. The AIC for the models under 
consideration is presented in the (Table |VlJ). It is clear that model which minimizing AIC 
is ACDM. Therefore it is no reason to introduce a model with bouncing terms and such 
model should be ruled out by Occam's razor. Because the extended bouncing dust model 
is statistically admissible from SNIa data it can be reconsidered only if the firm theoretical 
reason appears. Only this situation can justify consideration of the model with a small, but 
non-zero bouncing term. 

The existence of the oldest high-redshift extragalactic (OHReG) objects could be used as 
a test of the cosmological models (Table rVTTj) . The globular cluster analysis indicated that 
the age of the Universe is 13.4 Gyr 3^|. We demonstrate that the age of OHReG objects 
restricts the model parameter. As a criterion we take that the age of the Universe in a given 
redshift should be bigger than, at least equal, to the age of its oldest objects. With the 
assumption of H = 65 km/s MPc, the age of the universe on particular z for three class 
of models is calculated (Fig. |7J). This test admits the ACDM model with f2 m = 0.3. In 
this model, the age of the Universe is 14.496 Gyr. The BM model seems to be allowed from 
this test, however, that model predicts much longer age of the Universe (more than 20 Gyr) 
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TABLE V: Results of comparison of ACDM model with the extended bouncing models (m = 3) 
with fixed values n = 4 and n = 6. The result of statistical analysis for SNIa data from the best 
fit with minimum x 2 (denoted as BF) and from the likelihood method (denoted as L). The case of 
a fixed value of Q m fi is denoted as F. 



morlpl 




_7a n M Y 2 


J-.LJ.V_. Li J J. WW. 




u.ox 




1 _) 1 




0.34 — 


67 15 965 — 


L 




F0.30 — 


70 15 955 175 9 

v 7 . t V / J- " 7 • \J '.7 '. 7 i- 1 '.7 • '7 


BF 




F0.30 — 


70 15 945 — 


L 


ABCDM model 


n 31 n nn 


D 69 15 955 1 75 9 


BF 


with n = 4 


0.37 0.00 


0.65 15.965 


L 




F0.30 0.00 


0.70 15.955 175.9 


BF 




F0.30 0.00 


0.70 15.945 


L 


ABCDM model 


0.31 0.00 


0.69 15.955 175.9 


BF 


with n = 6 


0.34 0.00 


0.66 15.965 


L 




F0.30 0.00 


0.70 15.955 175.9 


BF 




F0.30 0.00 


0.70 15.945 


L 



then ACDM. In turn the BCDM model must be rejected because its age is 11.5 Gyr. 

IV. CMB PEAKS IN THE EXTENDED BOUNCING MODEL 

Acoustic oscillations in the primeval plasma during the last scattering give rise to the 
temperature map of cosmic microwave background (CMB). Peaks in the power spectrum 
correspond to maximum density of the wave. In the Legendre multipole space these peaks 
correspond to the angle subtended by the sound horizon at the last scattering. Further 
peaks answer to higher harmonics of the principal oscillations. 

The locations of these peaks depend on the variations in the model parameters. Therefore, 
they can be used to constrain the parameters of cosmological models. 
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FIG. 4: For the extended bouncing model with M. = 15.955) there are shown the confidence levels 
on the plane (fi raj o,n) minimized over parameter Q m ,o- 

The acoustic scale £4 which puts the locations of the peaks is denned as 

" 2 dcc dz' 



/c 



^0 H(z') , , 

Ea = roo ~ dz > (odj 
J 2. 



lz dec ^H(z>) 

where 



if (*) = ff ^ m ,o(l + Zf + a, (l + ZY - Q n , (l + *)» + Ov,0 ( 34 ) 

and c s is the speed of sound in the plasma given by 

c2 ^d M= ln Jfi (i + z) - ^ n n nfi (i + z) n - 3 

° s d PcS 3fi b)0 + 4Q 7;0 (1 + z)- nQ nfi (l + z) n ~ 3 ' 
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FIG. 5: Extended Bouncing model with M=15.955. The density distribution for £l n fi. Confidence 
level 68.3% and 95.4% are also marked on the figure. 

The properties of the bouncing term are unknown. In particular, we do not know 
whether it influences the sound velocity. But we assume that sound can propagate in it as 
well as in baryonic matter and photons. Let us note that with the lack of the bouncing term 
(i.e. fl n! o = 0) and/or when sound does not propagate in the bouncing fluid, we obtain the 
standard formula for Cg [4(1 ]. 

Knowing the acoustic scale we can determine the location of m-th peak 



4i ~ ^(m - (f) m ) 



(36) 



where <p m is the phase shift caused by the plasma driving effect. Assuming that f2 m) o = 0.3, 
on the surface of last scattering Zd ec it is given by 



0.267 



0.3 



o.i 



0.267 



1 Pr{Zda 

0.3 p m (z de 



0.1 



0.267 



1 f2 ri o(l + z de 
0~3 03 



o.i 



(37) 
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FIG. 6: Extended Bouncing model with M=15.955. The density distribution for n. Confidence 
level 68.3% and 95.4% are also marked on the figure. 

where fl hfi h 2 = 0.02, r(z doc ) = p r {z dcc ) / p m (z dec ) = Q T)Q (1 + z dec )/Q mfi is the ratio of the 
radiation to matter densities at the surface of the last scattering. 

The locations of the first two peaks are taken from the CMB temperature angular power 
spectrum j^l, 42|, while the location of the third peak is from the BOOMERANG measure- 
ments 0| . They values with uncertainties on the level la are the following 

4 = 220.11^, 4 = 546^°, 4 = 845+^. 

From the WMAP data, only the Hubble constant is H Q = 72 km/s MPc (or the parameter 
h = 0.72), the baryonic matter density i\,o = 0.024/i~ 2 , and the matter density fi m0 = 
0.14/i~ 2 [41] which give a good agreement with the observation of the position of the first 
peak. 

In the analysis of the constraints on the bouncing cosmological model parameters we fix 
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TABLE VI: The Akaike information criterion (AIC) for models under consideration: Einstein-de 
Sitter model (CDM), ACDM model (ACDM), bouncing model (BM), bouncing model with dust 
m = 3 (BCDM) and extended bouncing model with dust m = 3 (ABCDM). 



TTlodpl 


no. of parameters AIC 


CDM 


I 


325.5 


ACDM 

1 V V.. 1 '1V1 


2 


179.9 


BM 


4 


189.6 


BCDM 


q 
o 


Oil A 


BCDM with n = 4 


2 


230.6 


BCDM with n = 6 


2 


295.2 


ABCDM 


4 


183.9 


ABCDM with n = 4 


3 


181.9 


ABCDM with n = 6 


3 


181.9 


TABLE VII: The age 


of extragalactic 


objects. 



No object z age in Gys 



1 globular cluster 0. 13 — 15 

2 3C65 quasar 1.175 4.0 

3 LBDS 53W069 1.43 4.0 

4 LBDS 53W091 1.55 3.5 

the baryonic matter density f2b ; o = 0.05, the spectral index for initial density perturbations 
n = 1, and the radiation density parameter |40J 

^r,o = ^7,0 + ^u,o = 2.48/T 2 x 10" 5 + 1.7/T 2 x 10~ 5 = 4.18/T 2 x 10~ 5 (38) 

which is a sum of the photon f2 7i0 and neutrino Q Ut o densities. 

Assuming fi m ,o = 0.3 and h = 0.72 we obtain for the standard ACDM cosmological model 
the following positions of peaks 

i x = 220, £ 2 = 521, £ 3 = 821 

with the phase shift 4> m given by (pT7j). 
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FIG. 7: The age of the universe on particular z for three class of models: ACDM (middle curve), 
bouncing model BM (upper curve) and bouncing model with dust matter BCDM (lower curve). 
We marked age of 4 extragalactic object by star fTable lVTTj) . 

From the SNIa data analysis, it was found that the Hubble constant has a lower value. 
Assuming that H = 65 km/s MPc (or h = 0.65), we have f2 r ,o = 9.89 x 1CT 5 from eq. (|38|). 
For further calculations we take f2 r ,o = 0.0001. If we consider the standard ACDM model, 
with f2 nij0 = 0.3, Qb t Q = 0.05, the spectral index for the initial density perturbations n — 1, 
and h = 0.65, where sound can propagate in baryonic matter and photons, we obtain the 
following locations of first three peaks 

t x = 225, £ 2 = 535, £ 3 = 847. 

We note the difference between the observational and theoretical values in this case. We 
check whether the presence of the bouncing term f2 n moves the locations of the peaks. 
We do not know whether it influences the sound velocity, but we assume that sound can 
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TABLE VIII: Values of fi n o an d location of first three peaks. 



Model 


Hubble 


constant 






f\ 

A I 


x Z 




Extended Bouncing model n = 


= 4 


#0 


= 65 


km/s 


MPc 


3.0 x 10" 


-4 


217 


517 


816 






#0 


= 72 


km/s 


MPc 


2.86 x 10 


-4 


222 


526 


829 


Extended Bouncing model n = 


= 6 


H 


= 65 


km/s 


MPc 


1.4 x 10" 


10 


223 


530 


847 






H 


= 72 


km/s 


MPc 


1.3 x 10" 


10 


224 


530 


847 



propagate in it as well as in baryonic matter and photons. 

To obtain the bounce, n > 4 is necessary because the presence of the radiation term is 
required by the physics of primordial plasma in the recombination epoch. From the location 
of the first peak we obtain, the limit for f2 n term. In the case n = 5, with H = 72 km/s 
MPc, we obtain that Q ni o < 2 x 10~ n while for n = 6 we have that Q n> o < 2 x 10~ 17 . 

Please note that the special case n = 6 was analyzed for both values H = 65 km/s MPc 
and Hn = 72 km/ s MPc and the agreement with the observation of the location of the first 

n 

peak was obtained, also for the non-zero values of the parameter f2 nj o 30]. It means that 
both values of Hq are allowed from the CMB constraints for the case n = 6. However, this 
value is of order 1CT 10 . The results of calculations of the peak locations and the values of 
the parameter f2 n are presented in Table IVIIII In the special case n = 4, the bounce term 
scale like radiation, the existing of the bounce requires f2 n > ^r,o- in this case, we also 
obtain the agreement with the observation of the location of the first peak for the non-zero 
values of the parameter f2 n>0 (in order to 3 x 10~ 4 ). 

Finally, we analyze the models in which we assume that sound can propagate only in 
baryonic matter and photons. With H = 72 km/s MPc, in the case of n = 5, we obtain 
that fi n)Q < 2.2 x 10~ 8 while for n = 6we have that f2 n) o < 5 x 10" 14 . For the special case 
n = 4 we have that Q n £ = 2.3 x 10~ 4 . 

We have also calculated the age of the Universe in the ABCDM model. We find that 
the difference in the age of the Universe is smaller than 10 mln years for all values of Q n fi 
admissible by the CMB peaks location. So this model is admissible by the test of the age of 
the OHReG objects. 
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V. CONSTRAINT FROM THE BBN 



The observations of abundance of light elements is in good agreement with the prediction 
of the standard big-bang nucleosynthesis (BBN). It means that the BBN does not allow for 
any significant divergence from the standard expansion law, apart from the beginning of 
BBN to the present epoch. Therefore, any nonstandard terms included in the Friedmann 
equation should give only a negligible small change during the BBN epoch to render the 
nucleosynthesis process unchanged. 

It is crucial for the bouncing models to be consistent with BBN. These models have the 
nonstandard term Q n which scales like a~ n where n > 4. For example we analyze the cases 
n = 5 and n — 6. This additional term scales like (1 + z) n . It is clear that such a term gives 
rise to the accelerated Universe expansion if VL n p > 0. Going backwards in time, this term 
would become dominant at some redshift. If it happened before the BBN epoch then the 
radiation domination would never occur and the all BBN predictions would be lost. 

The domination of the bouncing term Q n should end before the BBN epoch starts and we 
assume that the BBN results are preserved in the bouncing models. In this way we obtain 
another constraint on the value of fi n - Let the model modification be negligible small 
during the BBN epoch and the nucleosynthesis process be unchanged. It means that the 
contribution of the bouncing term fi n cannot dominate over the radiation term Q r « 10~ 4 
before the BBN (z ~ 10 s ) 

|fi n ,oi(i + ^) n <a,o(i + ^) 4 . 

It means that |fi n ,o| < 10 -20 for the case n = 6 while |f) n o| < 10~ 12 for the case n = 5 
respectively. Of course, the case n = 4 is excluded because the existence of bouncing 
requires in this case \fl n ,o\ > ^r,o, while BBN constraints require \fl n ,o\ < ^r,o- Let us note 
that inequality Xf, < " 4 constrains the minimal size of the universe. The general 

conclusion from BBN constraints is that in the present epoch, the bouncing term, if it exists, 
is insignificant in comparison to the matter term. 

Table |IX] gives the value of bounce calculated for the best-fitted model parameters. Be- 
cause the bounce should take place before BBN epoch so ^bounce > ^bbn — 10 8 . Comparing 
with the -Zbounce presented in Table iDCl we obtain that only two classes of models ACDM and 
ABCDM are admissible. 



27 



TABLE IX: The value of ^bounce f° r the models under consideration: Einstein-de Sitter model 
(CDM), ACDM model (ACDM), bouncing model (BM), bouncing model with dust m = 3 (BCDM) 
and extended bouncing model with dust m = 3 (ABCDM). 



model 


ounce 


Vj\J Ivi 






AU1J1V1 






BM 




3.54 x 10 4 


BCDM 




2.98 


BCDM with n -- 


= 4 


3 


BCDM with n -- 


= 6 


4.05 x 10 4 


ABCDM 




oo 


ABCDM with n 


= 4 


oo 


ABCDM with n 


= 6 


oo 



VI. CONCLUSION 



In this paper we confront the bouncing model with astronomical observations. We use the 
constraints from SNIa data, CMB analysis, and BBN and the age of the oldest high-redshift 
objects. 

The standard bouncing model is excluded statistically at the 4cr level. If we take the 
extended bouncing model (with extra fi^o term) then we obtain, as the best-fit, that the 
parameter is equal to zero which means that the SNIa data do not support the exis- 
tence of the bouncing term in the model. We also demonstrate that BBN gives stringent 
constraints on the extra term and show that the bounce term is insignificant in the 
present epoch. 

It is interesting that such bouncing models with extra inflationary expansion are presently 
favored in the loop quantum approach (44 , The theory of loop quantum gravity 

predicts that there is no initial singularity because of the quantum effects in the Planck 
scale. It is due to the continuum break and granuality of spacetime. Therefore, we consider 
the model where we assume a small positive value of Q Ut o and estimate the rest of the 
parameters. This model is statistically admissible. However, when we compare this model 
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with the standard ACDM model applying the Akaike criterion, the latter is preferred. 

If the energy density is so large then quantum gravity corrections are important at both 
the big-bang and big-rip. It is interesting that the classical theory reveals its own bound- 
aries (i.e. classical singularities). The account of quantum effect avoids not only an initial 
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singularity but allows also to escape from a future singularity 

The avoidance of the initial singularity arises only on the quantum ground because the 
classical theory of gravity according to the Hawking- Penrose theorems states that these sin- 
gularities are essential if only some reasonable conditions on the matter content are fulfilled. 

If we assume the classical gravity is obvious during the whole evolution of the Universe 
than there is no reason to introduce the bouncing era. The ACDM model with the big-bang 
is a simpler model, while the bouncing model requires to the admittance of observationally 
unconfirmed assumptions. In this way Occam's razor methodology rules out the generalized 
bouncing model. The general conclusion is that the present astronomical data does not 
support the bouncing cosmology. 

We also adopt the methods of dynamic systems for investigating dynamics in the phase 
space. The advantages of these methods are that they offer the possibility of the investigation 
of all evolutional paths for all initial conditions. We show that the dynamics can be reduced 
to a two-dimensional Hamiltonian system. We also show structural instability of both the 
standard and generalized bouncing models. Let us note that the concordance ACDM models 
are structurally stable plj . The structural stability is a reasonable condition which should 
be satisfied by models of real physical processes. From the dynamic investigation we obtain 
that all models with the bounce are rather fragile. It means that any small perturbation of 
the right-hand sides of the dynamic equations of the model changes the topological structure 
of the phase space. The bouncing models in the space of all dynamic system on the plane 
form non-dense (zero measure) subset of this plane following the Peixoto theorem. Therefore, 
the bouncing models are untypical while ACDM models are generic from the point of view 
structural stability. 
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